as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="vessile")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ground")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ignited")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
View(mutated_data)
data<-read_csv("combind_data.csv")
pivoted_data<-data%>%pivot_longer(cols = c(dried,ignited,ground,raw,sived,vessile),
names_to = "Prosessing_Step",
values_to = "Conc.")%>%
#filter(Prosessing_Step%in%c("vessile","ground","ignited"))%>%
pivot_wider(id_cols = c(ID:M_Replicat,Prosessing_Step),
names_from = element,
values_from = Conc.)
library(tidyverse)
library(VCA)
data<-read_csv("combind_data.csv")
pivoted_data<-data%>%pivot_longer(cols = c(dried,ignited,ground,raw,sived,vessile),
names_to = "Prosessing_Step",
values_to = "Conc.")%>%
#filter(Prosessing_Step%in%c("vessile","ground","ignited"))%>%
pivot_wider(id_cols = c(ID:M_Replicat,Prosessing_Step),
names_from = element,
values_from = Conc.)
mutated_data<-pivoted_data%>%mutate(large_het=str_c(Depth,Sample_Replicat),
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:5,~as.factor(.)))
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ignited")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
mutated_data<-pivoted_data%>%mutate(large_het=str_c(Depth,Sample_Replicat))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
#relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:5,~as.factor(.)))
View(mutated_data)
mutated_data<-pivoted_data%>%mutate(large_het=str_c(Depth,Sample_Replicat))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
#relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:4,~as.factor(.)))
filterd_data<-mutated_data%>%
select(1:4,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ignited")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:4,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="sived")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:4,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="dried")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:4,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="vessile")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
mutated_data_2<-pivoted_data%>%
mutate(large_het=str_c(Depth,Sample_Replicat),
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:4,~as.factor(.)))%>%
filter(Prosessing_Step%in%!c("vessile","ground","ignited"))
mutated_data_2<-pivoted_data%>%
mutate(large_het=str_c(Depth,Sample_Replicat),
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:4,~as.factor(.)))%>%
filter(!Prosessing_Step%in%c("vessile","ground","ignited"))
View(mutated_data_2)
bind_cols(mutated_data_1,mutated_data_2)
mutated_data_1<-pivoted_data%>%
mutate(large_het=str_c(Depth,Sample_Replicat),
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:4,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
mutated_data_2<-pivoted_data%>%
mutate(large_het=str_c(Depth,Sample_Replicat))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(!Prosessing_Step%in%c("vessile","ground","ignited"))
mutated_data_2<-pivoted_data%>%
mutate(large_het=str_c(Depth,Sample_Replicat))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(!Prosessing_Step%in%c("vessile","ground","ignited"))
bind_cols(mutated_data_1,mutated_data_2)
mutated_data<-bind_cols(mutated_data_1,mutated_data_2)
View(mutated_data)
mutated_data<-bind_rows(mutated_data_1,mutated_data_2)
mutated_data_1<-pivoted_data%>%
mutate(large_het=str_c(Depth,Sample_Replicat),
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
mutated_data_2<-pivoted_data%>%
mutate(large_het=str_c(Depth,Sample_Replicat))%>%
select(!c(ID,Site,Date,Sive_Fraction,Depth,Sample_Replicat))%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:4,~as.factor(.)))%>%
filter(!Prosessing_Step%in%c("vessile","ground","ignited"))
mutated_data<-bind_rows(mutated_data_1,mutated_data_2)
View(mutated_data)
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="raw")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="sived")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="dried")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
View(mutated_data)
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/Instrument, method = "reml")
View(filterd_data)
View(mutated_data)
View(mutated_data)
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="vessile")%>%
as.data.frame()
View(filterd_data)
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/Instrument, method = "reml")
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ground")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ignited")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
levels(mutated_data$Prosessing_Step)
steps<-levels(mutated_data$Prosessing_Step)
steps<-levels(mutated_data$Prosessing_Step)
relm_outputs<-list()
for (prosessing in steps) {
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step==prosessing)%>%
as.data.frame()
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
relm_outputs[[prosessing]] <- model$aov.tab
}
steps<-levels(mutated_data$Prosessing_Step)
relm_outputs<-list()
for (prosessing in steps) {
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step==prosessing)%>%
as.data.frame()
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
relm_outputs[[prosessing]] <- model$aov.tab
}
View(relm_outputs)
steps<-levels(mutated_data$Prosessing_Step)
relm_outputs<-list()
for (prosessing in steps) {
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step==prosessing)%>%
as.data.frame()
if (prosessing%in%c("raw","sived","dried")) {
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
}else{
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
}
relm_outputs[[prosessing]] <- model$aov.tab
}
View(relm_outputs)
relm_outputs[[1]]
relm_outputs[[6]]
#creat vector and list for loop
steps<-levels(mutated_data$Prosessing_Step)
relm_outputs<-list()
for (prosessing in steps) {
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step==prosessing)%>%
as.data.frame()
if (prosessing%in%c("raw","sived","dried")) {
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het, method = "reml")
}else{
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/large_het/M_Replicat, method = "reml")
}
relm_outputs[[prosessing]] <- model
}
View(relm_outputs)
relm_outputs[[1]]$$aov.tab
relm_outputs[[1]]$aov.tab
relm_outputs[[1]]
relm_outputs[[1]]
#setup df for data with instrument replicats
mutated_data_1<-pivoted_data%>%
mutate(large_het=Sample_Replicat,
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(c(Prosessing_Step,large_het),.before = M_Replicat)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
View(mutated_data_1)
#setup df for data with instrument replicats
mutated_data_1<-pivoted_data%>%
mutate(large_het=Sample_Replicat,
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Sample_Prosessing,.before = Sample_Location)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
#setup df for data with instrument replicats
mutated_data_1<-pivoted_data%>%
mutate(large_het=Sample_Replicat,
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Prossesing_Step,.before = Sample_Location)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
#setup df for data with instrument replicats
mutated_data_1<-pivoted_data%>%
mutate(large_het=Sample_Replicat,
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Prossesing_Step,.before = Sample_Location)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
#setup df for data with instrument replicats
mutated_data_1<-pivoted_data%>%
mutate(large_het=Sample_Replicat,
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Prosessing_Step,.before = Sample_Location)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
View(mutated_data_1)
#setup df for data with instrument replicats
mutated_data_1<-pivoted_data%>%
mutate(large_het=Sample_Replicat,
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Prosessing_Step,.before = Sample_Location)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
#setup df for data without instrument replicats
mutated_data_2<-pivoted_data%>%
mutate(large_het=Sample_Replicat)%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Prosessing_Step,.before = Sample_Location)%>%
mutate(across(1:4,~as.factor(.)))%>%
filter(!Prosessing_Step%in%c("vessile","ground","ignited"))
#recombine data
mutated_data<-bind_rows(mutated_data_1,mutated_data_2)
#setup df for data with instrument replicats
mutated_data_1<-pivoted_data%>%
mutate(large_het=Sample_Replicat,
Instrument=if_else(M_Replicat %in% c(2, 4, 6), 2, 1),
M_Replicat=if_else(M_Replicat%in%c(1,2), 1,
if_else(M_Replicat%in%c(3,4), 2,
if_else(M_Replicat%in%c(5,6), 3,NA))))%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(Instrument,.before = P)%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Prosessing_Step,.before = Sample_Location)%>%
mutate(across(1:6,~as.factor(.)))%>%
filter(Prosessing_Step%in%c("vessile","ground","ignited"))
#setup df for data without instrument replicats
mutated_data_2<-pivoted_data%>%
mutate(large_het=Sample_Replicat)%>%
select(!c(ID,Site,Date,Sive_Fraction,Sample_Replicat))%>%
relocate(large_het,.before = M_Replicat)%>%
relocate(Prosessing_Step,.before = Sample_Location)%>%
mutate(across(1:5,~as.factor(.)))%>%
filter(!Prosessing_Step%in%c("vessile","ground","ignited"))
#recombine data
mutated_data<-bind_rows(mutated_data_1,mutated_data_2)
#creat vector and list for loop
steps<-levels(mutated_data$Prosessing_Step)
relms_outputs<-list()
for (prosessing in steps) {
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step==prosessing)%>%
as.data.frame()
if (prosessing%in%c("raw","sived","dried")) {
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het, method = "reml")
}else{
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
}
relms[[prosessing]] <- model
}
#creat vector and list for loop
steps<-levels(mutated_data$Prosessing_Step)
relms<-list()
for (prosessing in steps) {
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step==prosessing)%>%
as.data.frame()
if (prosessing%in%c("raw","sived","dried")) {
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het, method = "reml")
}else{
model<-fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
}
relms[[prosessing]] <- model
}
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="raw")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="sived")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="dried")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="vessile")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ground")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ignited")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
filterd_data<-mutated_data%>%
select(1:5,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ground")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
View(filterd_data)
filterd_data<-mutated_data%>%
select(1:6,Fe)%>%
filter(!is.na(Fe),
Prosessing_Step=="ground")%>%
as.data.frame()
fitVCA(Data = filterd_data, Fe ~ Sample_Location/Depth/large_het/M_Replicat, method = "reml")
View(filterd_data)
setwd("C:/Users/s299038/OneDrive - Cranfield University/1. Rapid Measurement Tools - Niall/Chapter 3 - XRF Analysis Accuracy on Slag/Data Analysis/Fig3 ANOVA")
#import data==========
setwd("C:/Users/s299038/OneDrive - Cranfield University/1. Rapid Measurement Tools - Niall/Chapter 3 - XRF Analysis Accuracy on Slag/Data Analysis/Fig3 ANOVA/source data")
# Get the list of CSV files in the working directory
csv_files <- list.files(pattern = "*.csv")
# Create an empty list to store the dataframes
data_tidys <- list()
i<-0
# Loop through each CSV file
for (file in csv_files) {
# Read the CSV file into a dataframe and pivot longer
df <- read_csv(file)%>%
filter(Sive_Fraction==1)%>%
mutate(across(Site:M_Replicat, ~ as.factor(.)),
across(P:"Pb +/-", ~ as.numeric(.) ),
Sive_Fraction=1,)%>%
relocate(ID)%>%
select(c(ID:Pb))%>%
select(!Vessile)%>%
pivot_longer(cols=P:Pb,
names_to = "element",
values_to = "Conc",)
list_name<-df%>% pluck("Prosesing_Step",1)%>%as.character()
df<-df%>%rename(!!list_name:=Conc)%>%
select(!Prosesing_Step)
# Add the dataframe to the list
data_tidys[[list_name]] <- df
#setup i
i<-i+1
#ifloop to join df together
if (i==1) {
joined_data<-data_tidys[[i]]
}else{
joined_data<-joined_data%>%
full_join(data_tidys[[i]],by=c("ID",
"Site",
"Date",
"Sample_Location",
"Depth",
"Sample_Replicat",
"Sive_Fraction",
"M_Replicat",
"element"))
}
joined_data%>%arrange(ID)
}
View(joined_data)
setwd("C:/Users/s299038/OneDrive - Cranfield University/1. Rapid Measurement Tools - Niall/Chapter 3 - XRF Analysis Accuracy on Slag/Data Analysis/Fig3 ANOVA")
